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Abstract 

A special normalization is proposed for strongly intensive quantities used in the study of event- 

o 

by-event fluctuations in high energy collisions. It ensures that these measures are dimensionless and 
yields a common scale required for a quantitative comparison of fluctuations of different, in general 
dimensional, extensive quantities. Namely, the properly normalized strongly intensive quantities as- 
sume the value one for fluctuations given by the model of independent particle production and zero 
in the absence of event-by-event fluctuations. 
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I. INTRODUCTION 



Intensive quantities are denned within the grand canonical ensemble of statistical mechanics. 
They depend on temperature and chemical potential(s), but they are independent of the system 
volume. Strongly intensive quantities [1] are, in addition, independent of volume fluctuations. 
They were suggested for and are used in studies of event-by-event fluctuations of hadron pro- 
duction in nucleus-nucleus collisions at high energies. This is because in these collisions, the 
volume of created states varies from collision to collision, and is difficult or even impossible to 
measure. 

Strongly intensive quantities are defined using two arbitrary, extensive state quantities A 
and B. Here, we call A and B extensive when the first moments of their distributions for the 
ensemble of possible states is proportional to volume. They are referred to as state quantities as 
they characterize the states of the considered system, e.g., final states (or equivalently events) 
of nucleus-nucleus collisions or micro-states of the grand canonical ensemble. For example, A 
and B may stand for multiplicities of pions and kaons in a particular state, respectively. 

The simplest family of strongly intensive quantities is given by the ratio of the first moments 
(i.e., average values) of A and B: 

R[A,B] = (1) 

where averaging (...) is performed over the ensemble of considered states. 

There are two families of strongly intensive quantities which depend on the second and first 
moments of A and B and thus allow the study of state-by-state fluctuations [lj . These are: 

A[A,B) = ^-[(B) u[A] - {A)u[B}] , (2) 
E[A,B) = ^\{B)u[A) + (A)u[B) - 2 ((AB) - (A){B)) 1 , (3) 



where 



»\a\ - {A2) - {A)2 u\m - {B2) ~ {B)2 u) 

"[A] = ^ , "[B] = ^ (4) 

are scaled variances of A and B. The normalization factors C*a and Cs are required to be 
proportional to the first moment of any extensive quantity. 



It is important to stress that A[A, B] and Y1[A, B] are independent of system size fluctuations 
not only for the grand canonical ensemble of states. They are also independent of the average 
number of sources and source number fluctuations in the model of independent particle sources, 
for example, in the wounded nucleon model [2]. 

Strongly intensive quantities for fluctuations have a long history. The first quantity of this 
type, introduced in 1992, was the so-called $ measure of fluctuations [3J. According to the 
current classification the $ measure belongs to the S family [TJ. It is defined as the difference 
of the quantity calculated for a studied ensemble (e.g., central Pb+Pb collisions) and its value 
obtained within an independent particle model (IPM) which preserves basic features of the 
ensemble. Thus, by construction, $ = if the studied ensemble satisfies the assumptions of the 
IPM. In general, $ is a dimensional quantity and it does not assume a characteristic value for 
the case of non-fluctuating A and B. The latter properties were clearly disturbing in numerous 
applications of $ when attempting to characterize fluctuations in experimental data [I] and 
models [5]. 

In this paper we propose a specific choice of the Ca and Cs normalization factors which 
makes the quantities A [A, B] and £[A, B] dimensionless and leads to A[A, B] = Y\A, B] = 1 
in the IPM. Moreover, from the definition of ALA, B] and E[A, B] it follows that A[A, B] = 
S[v4, B] = in the case of absence of fluctuations of A and B, i.e., for oo[A] = oj[B] = (AB) — 
(A)(B) = 0. Thus the proposed normalization of A[A, B] and £L4, B] leads to a common scale 
on which these fluctuation measures calculated for different state quantities A and B can be 
compared. 

The paper is organized as follows. In Section [IT] we introduce an independent particle model 
within which we calculate the A[A, B] and B] quantities. The calculation details are given 
in Appendix [Aj Appendix [B] gives explicit expressions of A [A, B] and T,[A, B] for three choices 
of the quantities A and B. Specific models which share the properties of the IPM are discussed 



in Section III Section IV presents the proposal for the normalization of A [A, B] and B], 
discusses the procedure to calculate them from a given ensemble of states and provides an 
illustration by numerical examples. A summary in Section M closes the article. 
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II. T\A,B\ AND A[A,B] IN AN INDEPENDENT PARTICLE MODEL 



The independent particle model assumes that: 

(i) the state quantities A and B (e.g., of a micro-state of the grand canonical ensemble, or 
of a final state of nucleus-nucleus collisions) can be expressed as 

A = on + a 2 + . . . + a N , B = ft + ft + . . . + ftv , (5) 

where <x,- and ft denote single particle contributions to A and B, respectively, and N is 
the number of particles; 

(ii) inter-particle correlations are absent, i.e. the probability of any multi-particle state is the 
product of probability distributions P(atj,{3j) of single-particle states, and these proba- 
bility distributions are the same for all j — 1, . . . , N and independent of iV : 

P N (a u Pi, a 2 , ft, ... , a N , ftv) = V(N) x P(a u ft) x P(a 2 , ft) x • • • x P(a N , ftv) , (6) 

where V(N) is an arbitrary multiplicity distribution of particles. 

It is easy to show (see Appendix [A]) that within the IPM the average values of the first and 
second moments of A and B are equal to: 

(A) = a(N), (A 2 ) = (N) + a 2 [(N 2 ) - (N)] , (7) 

(B) =0{N), (B 2 ) =W(N)+ f [(N 2 ) - (N)} , (8) 
(AB) = a/3 (N) + a ■ P [(N 2 ) - (N)] . (9) 

The values of (A) and (B) are proportional to the average number of particles (N) and, thus, 
to the average size of the system. These quantities are extensive. The quantities a, (5 and a 2 , 
(3 2 , a (3 are the first and second moments of the single-particle distribution P(a,(3). Within 
the IPM they are independent of (N) and play the role of intensive quantities. 

Using Eq. ^ the scaled variance uj[A] which describes the state-by-state fluctuations of A 
can be expressed as: 

(A 2 ) - (A) 2 ^ - a 2 (N 2 ) - (iV) 2 r n rArl , N 

u[A\ = K (A) = = + « (N) = + a "[N] , (10) 
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where u[a] is the scaled variance of the single-particle quantity a, and u[N] is the scaled 
variance of N. A similar expression follows from Eq. ^ for the scaled variance u[B). The 
scaled variances w[A] and u[B] depend on the fluctuations of the particle number via w[iV]. 
Therefore, u)[A] and oo[B] are not strongly intensive quantities. 

From Eqs. ([7j|9]) one obtains expressions for A [A, B] and S[A, B\ 1 namely: 

(N) 



A[A,B] = 
E[A,B] = 
Thus, the requirement that 



(TV) 



(3 u[a] — a u[(3] 
P w[a] + a uj[(3] - 2 ( ~a]3 - a ■ p ) 



(11) 
(12) 



A[A, B] = Z{A,B] 



(13) 



within the IPM leads to: 



(N) pu[a] - au[p] 



(N) P<u[a] + auj[p] - 2{ap-a-p) 



(14) 
(15) 



Two comments are in order here. First, Eqs. (|7jj9]) have the same structure as Eqs. (2-4) 
of Ref. pQ obtained within the model of independent sources. The only difference is that 
the number of sources N$ in the model of independent sources is replaced by the number of 
particles iV in the IPM. Each source can produce many particles, and the number of these 
particles varies from source to source and from event to event. Besides, the physical quantities 
for particles emitted from the same source may be correlated. Therefore, in general, the model 
of independent sources does not satisfy the assumptions of the IPM. Nevertheless, the formal 
similarity between the two models can be exploited and gives the following rule of one to one 
correspondence: all results for the IPM can be found from the expressions obtained within the 
model of independent sources, assuming artificially that each source always produces exactly 
one particle. Second, only the first and second moments of two extensive quantities A and B 
are required in order to define the strongly intensive quantities A and E. However, in order 
to calculate the proposed normalization factors Cs and Ca additional information is needed, 
namely the first and second moments of single-particle contributions to A and B as well as the 
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mean number of particles. Note that in special cases the factors C^, and C*a may assume the 
value zero and thus the proposed normalization is not possible. 



Explicit expressions for Eqs. (14, 15) for three choices of A and B are given in Appendix [Bj 
The first two cases correspond to the study of "transverse momentum" and "chemical" fluctu- 
ations. The third choice is the most general. 

III. EXAMPLES OF INDEPENDENT PARTICLE MODELS 

In this section two specific models which satisfy the IPM assumptions, i.e., Eqs. ([5j [6]), are 
presented and discussed. 

A. Grand Canonical Ensemble 

The most popular model which satisfies the IPM assumptions is the ideal Boltzmann multi- 
component gas in the grand canonical ensemble formulation. Here we refer to it as the IB-GCE. 
In the IB-GCE the probability of any microscopic state is equal to the product of probabilities 
of single-particle states. These probabilities are independent of particle multiplicity. Thus, the 
IB-GCE satisfies the assumption ^ of the IPM. 

The IB-GCE predicts a specific form of the multiplicity distribution V(N), namely the 
Poisson distribution and thus, u[N] = 1. Moreover, it also predicts the specific form of the 
single-particle probability in momentum space, namely the Boltzmann distribution: 



Mp) = Cexp|- v/p2 T +m j , (16) 

where p and m are particle momentum and mass, respectively, T is the system temperature, 
and C — \J d 3 p fs(p)] 1 is the normalization constant. 

Note that by introducing quantum statistics one destroys the correspondence between the 
GCE and the IPM. This is because of (anti-) correlation between particles in the same quantum 
state for the (Fermi) Bose ideal gas. Moreover, correlations between particles are introduced if 
instead of resonances their decay products are considered. Note that it is necessary to include 
the strong decays of resonances in order to compare the GCE predictions to experimental 
results. 
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The correspondence between the IB-GCE and the IPM remains valid even if the volume 
varies from micro-state to micro-state 1 but local properties of the system, i.e., temperature and 
chemical potentials are independent of the system volume. Let volume fluctuations be given by 
the probability density function F(V). The averaging over all micro-states includes the aver- 
aging over the micro-states with fixed volume and the averaging over the volume fluctuations. 
The volume fluctuations broaden the V(N) distribution and increase its scaled variance: 



where V k = J dV F(V) V k for k = 1,2. The first term on the right hand side of Eq. (17) 
corresponds to the particle number fluctuations in the IB-GCE at a fixed volume V (i.e., this 
is the scaled variance of the Poisson distribution), and the second term is the contribution due 
to the volume fluctuations. Equation (J6| remains valid in this example, therefore, the IB-GCE 
with arbitrary volume fluctuations satisfies the IPM assumptions. 



B. Mixed Event Model 



The Mixed Event Model is defined by the Monte Carlo procedure frequently used by experi- 
mentalists in order to create a sample of artificial events in which correlations and fluctuations 
present in the original ensemble of events are partly removed. Then the original and mixed 
events are analyzed in the same way and the corresponding results are compared in order to 
extract the magnitude of a signal of interest, which by construction should be present in the 
original events and absent in the mixed events. The mixed event procedure is in particular 
popular in studies of resonance production, particle correlations due to quantum statistics and 
event-by-event fluctuations, see for examples Ref. [7]. 

There are many variations of the Mixed Event Model. Here we describe the one which in 
the limit of an infinite number of the original and mixed events gives results identical to the 
IPM. 

The procedure to create a mixed event which corresponds to the given ensemble of original 
events consists of two steps, namely: 

1 The statistical ensembles with volume fluctuations were discussed in Ref. [B] 
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(i) a mixed event multiplicity, N, is drawn from the set of multiplicities of all original events; 

(ii) N particles for the mixed event are drawn randomly with replacement from the set of all 
particles from all original events. 

Then the steps one and two are repeated to create the next mixed event and the procedure is 
stopped when the desired number of mixed events is reached. In the limit of an infinite number 
of original events, the probability to have two particles from the same original event in a single 
mixed event is zero and thus particles in the mixed events are uncorrelated. Therefore, in this 
limit the mixed event model satisfies the IPM assumptions. Note, for an infinite number of 
mixed events, the first moments of all extensive quantities and all single-particle distributions 
of the original and mixed events are identical. 



IV. NORMALIZATION AND DETERMINATION OF A[A,B] AND E[A,B] 

The strongly intensive quantities A[A, B] and E[A, B] were introduced for the study of state- 
by-state fluctuations of any extensive quantities A and B in a given ensemble of states. For 
example, states may refer to data for nucleus-nucleus collisions recorded by an experiment 
or generated within a Monte-Carlo model. In this section, we first explicitly introduce our 
special normalization of ALA, B] and E[A, B] for the ensemble of states. Then, we discuss the 
procedure of their determination. 

Let the ensemble of states Q and the extensive state quantities A and B be given. We propose 
to define the normalization factors Cs and Ca in Eqs. (|2j [3]) such that A [A, B] = E[A, B] = 1 
in the IPM with the multiplicity distribution V(N) and the single-particle distribution P(a,fl) 
identical to those of the ensemble Q. The IPM which corresponds to the ensemble Q will be 
denoted as the IPM-f2. The normalization factors Ca and Cs calculated in the IPM-f2 are then 



given by Eqs. (14, 15), where all entering quantities should be calculated from the ensemble f2. 



The procedure of calculating A [A, B] and E[A, B] given by Eqs. pi pi) with the normalization 



factors defined by Eqs. (14, 15) consists of the following steps: 



(i) calculate the f2-ensemble state averages of the first and second moments of extensive 
quantities A and B; 



8 



.— . 4.5f 
3 3.5 - 



UrQMD 3.3: 
□ P+P 
■ Xe+La (7%) 



3^ 
2.5 r 

2^ 
1.5-" 

0.5^ 

0^ 



-Q- 



•Q— H-- 



-a 



i ... i ... i ... i ... i 



6 8 10 12 14 16 18 20 



s NN [GeV] 



> 

(D 



Q. 

3 1 








UrQMD 3.3: 






(b) 


2.5 


r □ 


P+P 








2 


■ 


Xe+La (7%) 








1.5 






.M-"" 






1 

0.5 


: w 

:B- 


M—Wl"' 

-Q---0- 











" i ■ 


, I , , , I , , 


. I . . . I . 







8 10 12 14 16 18 20 



'NN 



[GeV] 



.— . 2 
"3: 1.8 



1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




r UrQMD 3.3: 






(c) 


: ■ P+P 




* 




L * Xe+La (7%) 


* 






















- 1 , , , i , , , i , 









8 10 12 14 16 18 20 



'NN 



[GeV] 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 




UrQMD 3.3: 
■ P+P 

★ Xe+La (7%) 




'NN 



[GeV] 



Figure 1: Fluctuation measures calculated within the UrQMD model for negatively charged hadrons 
produced in inelastic p+p interactions and the 7% most central Xe+La collisions as functions of 
collision energy in the CERN SPS energy range. The top plots show intensive measures of fluctuations, 
namely the scaled variances of (a) the negatively charged hadron multiplicity, uj[H~], and (b) the sum 
of magnitudes of their transverse momenta, cjfPy]. The bottom plots show the corresponding strongly 
intensive measures (c) A[Py, H~] and (d) Statistical uncertainties are smaller than the 

symbol size and were calculated using the subsample method. 
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(ii) calculate the first and second moments of single-particle quantities, a and /3, as well as 



the average number of particles (N) entering Eqs. (14, 15); the avareging is performed 
over the f2-ensemble states and particles; 

(iii) calculate A[A,B] and E[A, B] according to Eqs. @ and Q. 

The proposed procedure is illustrated by numerical results obtained within the UrQMD 
model [8] . Figure [l] shows the collision energy dependence of different fluctuation measures 
discussed in this paper in the CERN SPS energy range. In this example, A = Pt is the 
number of negatively charged hadrons and B = H is their total transverse momentum (see 



Appendix B 1 ). The UrQMD simulations were performed for inelastic p+p interactions and for 
the 7% most central Xe+La collisions. This choice of reactions is motivated by the experimen- 
tal program of the NA61/SHINE Collaboration [9j at the CERN SPS. NA61/SHINE already 
reported the first results on event-by-event fluctuations in p+p interactions [10], and results 
for nucleus-nucleus (Be+Be, Ar+Ca, Xe+La) collisions will become available within the next 
couple of years. A comparison between experimental data and models is beyond the scope 
of this paper. The top plots show intensive fluctuation measures, namely the scaled variance 
of the negatively charged particle multiplicity distribution, u[H~], and of the distribution of 
the sum of the magnitudes of their transverse momenta, u;[Pr]- The bottom plots show the 
corresponding strongly intensive measures A[Pt,H~] and £[Pr,P~] normalized as proposed 



in this paper according to Eqs. (14, 15) with their explicit form given in Eq. (B5). 

The scaled variance of H~ and Pt is significantly larger in central Xe+La collisions than 
in p+p interactions. To a large extent this is due to fluctuations of the number of nucleons 
which interacted (wounded nucleons), see Ref. [H] for a detailed discussion of this issue. The 
advantages of the A[Pt, H~] and E[Pr, H~] quantities are obvious from the results presented in 
the bottom plots. First, they are not directly sensitive to fluctuations of the collision geometry 
(the number of wounded nucleons) in contrast to the scaled variance. Thus, the remaining 
small differences between results for central Xe+La collisions and p+p interactions are entirely 
due to deviations of the UrQMD model from the independent source model. Second, they 
are dimensionless and expressed in units common for all energies and reactions as well as for 
different choices of state quantities A and B. Due to the particular normalization, proposed 
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in this article, they assume the value one for the independent particle model and zero in the 
absence of event-by-event fluctuations. 

V. SUMMARY 

Strongly intensive quantities A[A, B] and S[A, B] are fluctuation measures which are in- 
dependent of the system volume and its fluctuations within the grand canonical ensemble of 
statistical mechanics. Moreover, they are independent of the number of wounded nucleons and 
its fluctuations within the Wounded Nucleon Model. Strongly intensive quantities are expected 
to be useful in studies of fluctuations in hadron production in nucleus-nucleus collisions at high 
energies. In this paper a special normalization of strongly intensive quantities is proposed. It 
ensures that they are dimensionless and yields a common scale enabling a quantitative com- 
parison of fluctuations of different extensive state quantities. With the proposed normalization 
the discussed strongly intensive quantities assume the value one for fluctuations given by the 
independent particle model and zero in the absence of state-by-state fluctuations. 

The paper includes details of calculations, explicit formulas for "transverse momentum" and 
"chemical" fluctuations as well as for the most general case of fluctuations of two extensive 
motional quantities for partly overlapping sets of particles. Moreover numerical examples are 
given using final states of high energy collisions generated by the UrQMD model. 
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Appendix A: Calculation details of the IPM 



In this Appendix details of the derivation of Eqs. (|7jj9j) within the Independent Particle 
Model defined by Eqs. ([5j [6]) are given. 

The functions entering Eq. ^ satisfy the normalization conditions: 



JV 



The average values of ay and a 2 are: 

= (a 2 ) = • • • = 
(a?> = <c^> = ...= 

Similarly, 

</5i> = </5 2 > = - - - = 
(/?i 2 ) = (/3 2 2 ) = • • • = 

At i 7^ j one finds: 



(ajv) 



dad/3 P(a,(3) = 1 . 

y dad/3 a P(a, (3) = a, 
dad (3 a 2 P '(a, (3) = a 2 



(/?jv> 



dad(3 f3 P(a,f3) = (3 , 
dad/3 /3 2 P(a,/3) = ~W ■ 



(aiaj) = {a i ){a j ) = a- a = a 2 , (frfy) = (/%)(/%) = (3 ■ (3 = (3 



The state averages of A and B are equal to: 



JV 



JV 



(A) = <J>;> = ^^E^i) = EW 5,JV = 

j=l JV j=l JV 

JV JV 

(5) = (I>> = EwE<ft) = Ew^ ,JV = 



j=l V j=l JV 

For the second moments of A and B one obtains: 

JV 

(A 2 ) = ( («! + a 2 + • • • + aw) 2 ) = ( a ] + ^E aiQt i ^ 

JV 



Ew 



JV 



5E<«i> + $E ^ iQ! j) 

j=l i<-i^i< ^ 



a 2 (A) + a 2 [(A 2 ) - (A)] 



(Al) 



(A2) 
(A3) 



(A4) 
(A5) 



(A6) 



(A7) 
(A8) 



(A9) 
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and similarly 

(B 2 ) = W(N) + f [(N 2 )-(N)] . (A10) 

Finally, for (AB) one finds: 

(AB) = ( (ax + a 2 + . . . + a N ) x (ft + /3 2 + . . . + /3 N ) ) 



A' 



' TV 



a 



/3<JV) + a-P [(N 2 )-(N)} . (All) 



Appendix B: Examples for three choices of A and B 

In this Appendix explicit expressions for A [A, B] and T,[A, B] and their normalization factors 
C*s and Ca calculated within the IPM are given for two popular choices for the extensive state 
quantities A and B which correspond to the study of "transverse momentum" and "chemical" 
fluctuations. Finally the most general case is considered, which corresponds to the selection of 
two extensive motional quantities for partly overlapping sets of particles. 

1. "Transverse momentum" fluctuations 

The first [3] and the most popular [U |5] application of the $ measure was the study of 
transverse momentum fluctuations. In the formalism, introduced in Ref. pQ, this corresponds 
to the following choice of the extensive state quantities A and B: 

A = P T = pfi+p® + . . . +p { T N) , (Bl) 
B = N = w w + w {2) + ...+ w w , (B2) 

where p^) is the absolute value of the transverse momentum of j th particle 2 and is the 
particle identity [12J which equals to one for all particles: w^' = 1. 
Thus, for the single-particle quantities: 

a = pt , (3 = w = 1 , (B3) 



Similarly one can consider sums of any other motional variables, e.g., particle energies, rapidities, etc. 
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one gets: 

a = pr , a 2 = i (3 = (3 2 = w = w 2 = l, a(3 = fFf , (B4) 
where p~r and p T are the average values of pt and p\ calculated from the properly normalized 



single-particle transverse momentum distribution 3 . Consequently, Eqs. (14, 15) give 



2 2 

Ca = C s = (iV) • Pt L PT = (N) ■ co[p T } . (B5) 

Pt 

As was already mentioned in Section [TTJ only the first and second moments of two extensive 
quantities Pt and N are required to calculate the strongly intensive measures A[Pr,N] and 
H[Pt,N]. However, in order to calculate the proposed normalization factors Ca and Cs ad- 
ditional information may be necessary. In the considered example, one also needs the second 



moment of the single-particle p^-distribution p T to calculate (B5). 

Let us recall here that Y,(Pt,N) is directly related to the $ PT measure of transverse mo- 
mentum fluctuations, for the explicit expression see Ref. pQ. The only difference is in the scale 
used to quantify fluctuations measured by both quantities. Namley, the $ PT measure is defined 
as the difference of the event quantity calculated for the studied ensemble (e.g., central Pb+Pb 
collisions) and its value obtained within the independent particle model. Consequently, $ PT = 
if the studied ensemble satisfies the assumptions of the IPM. Moreover, $ PT is a dimensional 
quantity and does not assume a characteristic value for the case of non-fluctuating A and B. 
These undesired properties of $ PT are removed when fluctuations are measured using S(Py, N) 
normalized as proposed in this article. 



2. "Chemical" fluctuations 



In the jargon of high energy nuclear physics "chemical" fluctuations refer to fluctuations 
of particle-type composition of the system. In order to be specific let us consider relative 
fluctuations of the number of charged pions 7r = 7r + + 7r~ and kaons K = K + + K~: 

A = K = w V +w W + ... +w ^\ B = tt = w W +w f) + ... +w W , (B6) 
3 In high energy physics single-particle distributions are called inclusive distributions. 
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where and w K are the pion and kaon identities of the j th particle 4 . Particle identities were 
introduced first in Ref. [T2] and used in the study of "chemical" fluctuations in terms of the $ 
measure jU EJ [12] . 

In this example one defines the kaon w K and pion identities as: Wf) = 1 and = 
if the j th particle is a kaon, and w K = and w-i? = 1 if the j th particle is a pion. 
For the single-particle quantities, 



a 



w K 



(B7) 



one obtains: 



Wk = w 



K 



(K)_ 
(N) 



= k 



w n = wi 



(N) 



1 — k , wk w-k = 



where N = K + n. Then from Eq. (B8) follows: 

-2 



U)[W K ] = 



w K - wk 



WkW^ 



w k 



= 1 — k , uj [w n ] 
k- (1 - k) . 



wl - wJ 2 



w n 



(B8) 



(B9) 
(BIO) 



Therefore, Eqs. (14, 15) give: 



C A = (N)-{l-2k) = (tt) - (K) 
C s = (N) = (vr) + (K) . 



(BH) 
(B12) 



As seen from Eqs. (Bll, B12) the normalization factors C*a and Cs depend only on the first 



moments of the extensive state quantities K and 7r. 

However, in general more information is needed to calculate C*a and As an illustration let 
us consider partly overlapping sets of particles, e.g., the number of charged kaons K = K + + K~ 
and all negatively charged particles H . The extensive state quantities A and B are: 



A = K 
B = H' 



w 



(i) 

K 



(2) , 
W K ' + . . . + W 



(N) 
K i 



(1) , (2) , , (AO 

W_ ' + W_ ' + . . . + W_ ' 



(B13) 
(B14) 



4 Similarly one can consider sums of any other particle identities, e.g., negatively charged particle, baryons, 
etc. 
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where and w_ are the kaon and negatively charged hadron identities of the j particle. 



,00 



th 



The kaon and negatively charged hadron wf. } identities are defined as: w K ! = 1 and 



U) 



,C0 



w 



U) 



if the j th particle is a K + , = 1 and w ( l' = 1 if the j th particle is a if , and 



U) 



(j) 



and w^P = 1 if the j th particle is a negative hadron but not a K 



:th 



K 



For the single-particle quantities, 



a 



w K 



one obtains: 



w K 



W- 



w 



(K+) + (K-) 



w 



(N) 



(N) 



= h. 



k+ + k_ 



= k , 

(K~) 
(N) 



where N = K + + H . Then from Eqs. (B16, B17) it follows: 

-2 



U\W K \ = 



W 2 K - w k~ 

wk ■ wZ 



fc_ 



1 - k , 
- k-h. 



wz 



U\W-\ 



1 - h- 



W- 



(B15) 

(B16) 
(B17) 

(B18) 
(B19) 



WkW- 

Therefore, for Eqs. ( T4p5[ ) one finds: 

C A = (N)-(h--k) = (H~) - (K) , (B20) 
C s = (N) ■ {h- + k-2kJ) = (H~) + (K) - 2(K~) . (B21) 

Thus in this case the normalization factors depend on (K) and (H~), and in addition on (K~). 



3. The most general case 

The most general case, which up to now was not considered in the literature, concerns 
relative fluctuations of two motional extensive quantities, e.g., energy of charged kaons Ek and 
transverse momentum of all negatively charged hadrons . These two sets of particles are 
partly overlapping. This example corresponds to the following choice of the extensive state 
quantities A and B: 

A = Ek = w$eM + w^e® + ... + w™*™ , (B22) 
B = P~ = w^pW + w & p f > + . . . + wWjJW , (B23) 
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where and w_ are the kaon and negatively charged hadron identities of j particle, and 



,(i) 



tii 



and p[^ are its energy and transverse momentum. Note that for e = pt = 1 Eqs. (B22) and 



(B23) are reduced to Eqs. (B13 B14), respectively. 



For the single-particle quantities, 



w K e 



P 



W-Pt 



one obtains: 

a = k • e , a 2 = k ■ e 2 , = h- - pi , (3 2 = h- ■ p 2 , a(3 = k- ■ ~epi 
where (n = 1, 2) 

^ = J dee n f K (e) , = J dp t p n t f(p t ) , 

Wt = J d3 P \/p 2 + m KPtfK-(p) ■ 



(B24) 

(B25) 

(B26) 
(B27) 



In order to calculate the averages (B26, B27) one needs to know the single-particle e-distribution 



for kaons /ft-(e), the ^-distribution for negatively charged hadrons f(pt), and the p-distribution 



for K . Then from Eq. (B25) follows: 

u[a] = u[e] + (l-k)-e, co[f3] = co[p t ] + (l-h-)-pi 
a(3 — a ■ (3 = k- ■ ~ep~ t — kh^e - pi . 

Finally one finds for the normalisation factors: 



C A = (Pt) ■ 


w[e] + p t 


- (E K ) ■ 


u[p t ] + e 


> 


= (Pt) ■ 


w[e] + Wt 


+ (E K ) ■ 


u)\pt] + e 





For the special case e = p t = 1 one gets: 

e = e 2 = p~ t = p\ = epl = 1 , 

leading to: 

u>[e] = u\pt] = , (E K ) -+ (K) , (P T ) -+ (H~) , 



(B28) 
(B29) 

(B30) 
(B31) 

(B32) 

(B33) 



and Eqs. (B30 B31) reduce to Eqs. (B20, B21), respectively. 
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